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FOREWORD 


The  work  covered  in  this  report  was  performed  in  the  FBM 
Geoballistics  Division,  with  the  exception  of  Appendix  B which 
was  contributed  by  Christopher  Gracey  of  the  Aeromechanics  Branch 
of  the  Exterior  Ballistics  Division.  Both  divisions  are  in  the 
Strategic  Systems  Department. 

This  report  was  reviewed  by  Carlton  W.  Duke,  Jr.,  Head,  FBM  Geo- 
ballistics  Division. 
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METRIC  CONVERSION  TABLE 


To  convert  from 


To Multiply  by 


feet  (ft) 

feet/second  (ft/sec) 
degrees  (angle) 


meters  (m) 

meters /second  (m/sec) 
radians  (rad) 


0.30480 

0.30480 

0.01745 
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INTRODUCTION 


In  this  report  equations  are  derived  for  guiding  a variable  lift 
reentry  vehicle  to  impact  with  an  earth  fixed  target.  The  concept 
of  a dive  line  located  at  the  target  is  used  to  meet  preselected  ter- 
minal conditions  on  impact  flight  path  angle  and  approach  azimuth. 

With  a simplifying  assumption  the  resulting  equations  can  be  put  in 
a cross  product  form. 

Preliminary  simulation  results  are  given  but  the  report  is  circu- 
lated primarily  to  stimulate  discussion  of  the  ideas  involved  and  ! 

demonstrate  feasibility  of  the  equations  rather  than  present  an  in- 
depth  design  study. 

The  next  section,  the  primary  one,  contains  a derivation  of  the 
guidance  equations.  This  is  followed  by  a short  section  of  Simula-  . 

tion  results  and  discussion  and  finally  two  appendices.  The  first  > 

describes  the  coordinate  frames  and  notation  used.  The  second,  by  i 

Christoper  Gracey,  contains  a derivation  of  a simplified  cross  pro- 
duct form  of  the  guidance  law. 
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DERIVATION  OF  THE  GUIDANCE  EQUATIONS 


The  guidance  law  derived  in  this  section  is  a modification  of  a 
technique  originated  fay  the  author  in  the  early  sixties  for  ascent 
guidance.  The  problem  to  be  solved  is  as  follows:  given  a dive  line 
at  the  target,  control  the  vehicle  so  that  at  target  impact  the  vehi- 
cle's velocity  vector  is  parallel  to  the  dive  line.  The  solution  is 
to  drive  the  horizontal  plane  and  vertical  plane  components  of  the 
line  of  sight  which  are  orthogonal  to  the  dive  line,  and  their  time 
rates  of  change,  to  zero  by  impact  time. 

The  procedure  to  accomplish  this  is  simple.  Two  components  of 
the  commanded  lift  vector  are  prescribed  as  functions  of  the  esti- 
mated time  to  go.  The  functions  have  a total  of  four  undetermined 
parameters.  These  are  chosen  so  that  when  time  to  go  is  zero  the 
line  of  sight  components  of  interest  and  their  derivatives  are  zero. 
The  third  component  is  determined  by  imposing  the  physical  constraint 
that  the  lift  vector  is  normal  to  the  vehicle's  relative  velocity 
vector. 

Usually,  derivatives  of  the  vehicle  motion  with  respect  to  in- 
ertial space  are  the  ones  available.  In  deriving  the  guidance  law, 
however,  derivatives  are  needed  with  respect  to  an  earth  fixed  frame. 
These  are  related  to  inertial  derivatives  by  the  Coriolis  formula 
which  is  used  without  comment  where  needed.  In  the  following,  a dot 
over  a vector  indicates  its  time  derivative  with  respect  to  inertial 
space,  while  a prime  indicates  differentiation  relative  to  the  earth. 
For  scalars  the  distinction  does  not  exist  and  either  symbol  may  be 
used. 

The  first  step  in  the  derivation  is  to  compute  components  of  the 
line  of  sight  and  their  derivatives  relative  to  the  dive/tum  frame. 
Appendix  A contains  definitions  of  the  quantities  used  in  the  follow- 
ing. 


The  line  of  sight  vector  r which  is  defined  as 
r - RT  - R 

can  be  expressed  in  component  form  as 


1 


from  which  it  follows 


r'  * rDUD  + rPUP  + rNUN 


r"  * r"U  + r"U  + r"U 
r rD  u rPUP  rNUN 


From  the  definitions  of  E,  H,  F and  W given  in  the  appendix  it 
follows  that 

E = r • U„ 


H - E'  - x'  ‘ Ur 


F - r • Uv 


W *>  F'  = x'  * U* 


The  derivative  x'  is  computed  by 


x'  * a)  x R - R 


The  second  derivative  r"  will  be  needed  also  and  is 


r"  ■ 2id  x R - 10  x (u>xR)  - R 


At  this  point  it  should  be  noted  that  r"  is  used  only  in  the 
closed-loop  guidance  equations;  so  approximations  are  in  order.  The 
general  philosophy  in  this  report  regarding  approximations  holds  that 
they  are  acceptable,  provided  the  errors  engendered  are  nowhere  too  bad, 
and  that  they  get  smaller  (or  remain  negligible)  as  terminal  conditions 
are  met. 

With  this  in  mind  consider  the  magnitudes  of  the  first  two  terms 
in  the  equation  for 


.1  _5  4 2 

1 2u>  x R I — 2 x (7x10  ) x (1x10  ) ■ 1.4  ft/sec 


3 


J 


I 


|w  x (wxR) \~-  (7x10  x 2x10^  = 0.1  ft/sec^ 


i 

Both  of  these  are  small  compared  with  the  available  acceleration  due 
to  lift;  so,  in  all  that  follows,  the  approximation  is  made. 


r 


II 
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Consistent  with  this  approximation, 


E"  = r"  • UD  = - R • U,. 


F"  =*  r"  * Ux,  - - R ’ U. 


The  next  step,  the  pivotal  one,  is  to  require 


E"  « A + BTg 
F"  = C + DTg 


where  TG  is  an  estimate  of  time  to  go. 
Three  problems  remain  to  be  solved: 


(a)  Determine  the  parameters  A,  B,  C and  D by  requiring 
satisfaction  of  the  desired  terminal  conditions. 

(b)  Find  an  approximation  for  time  to  go. 

(c)  Relate  components  of  the  lift  vector — the  only  quantity 
available  for  control — to  the  parameter  values  and  Tg. 

If  all  the  approximations  made  were  exact.  A,  B,  C and  D would  not  vary 
with  time.  In  practice  they  do  vary  under  continual  recomputation,  but 
this  is  ignored  when  computing  their  values  at  any  time  step. 

In  the  dive/turn  frame,  the  lift  vector  is  expressed  as 


LU 
D D 


+ LU  +LU 
P P NN 
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The  three  components  are  determined  by  imposing  the  constraints, 

(1)  E and  E'  must  approach  zero  as  TG  goes  to  zero 

(2)  F and  F'  must  approach  zero  as  T goes  to  zero 

G 

(3)  The  lift  vector  is  perpendicular  to  the  relative 
velocity  vector. 

To  impose  condition  (1) , the  equation  for  E"  is  integrated  twice 
with  respect  to  TG,  noting  the  sign  change  and  taking  into  account 
current  values  of  E and  E'.  This  procedure  yields  two  equations  for 
A and  B,  viz. 

2A  + BTG  = - 2H/Tg 
2 

3A  + BT_  - GE/T 

Vj  G 

which  have  solutions. 


In  exactly  the  same  way,  condition  (2)  implies 


Before  implications  of  condition  (3)  are  considered,  it  is  perhaps 
as  well  to  digress  for  a moment  to  find  an  expression  for  Tq  and 
show  that  the  constraints  imposed  so  far  have  determined  Lp  and  Lg. 
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First,  an  estimate  for  time  to  go  is  obtained  simply  by  dividing 
the  magnitude  of  the  line  of  sight  vector  r"  by  the  negative  of  its 
time  rate  of  change.  This  yields 


T 


G 


r • r 
r'  • T 


or,  equivalently. 


where  V is  the  relative  velocity  vector 


Next,  to  show  that  conditions  (1)  and  (2)  determine  Lp  and  LN, 
consider  the  total  acceleration  of  the  vehicle  with  respect  to 
inertial  space.  It  is  just  the  sum  of  the  accelerations  due  to  lift, 
drag  and  gravitation. 

Since  drag  changes  the  direction  of  the  relative  velocity  little 
compared  with  lift,  it  will  be  neglected  in  computing  the  commanded 
! lift.  With  this  comment,  r1'  becomes 

I 

| r"  - - (L  + G) 

where  G is  the  gravitation  vector.  It  follows  immediately  that 


Lp  - - (A  + BTg) 

Ln  - - (C  + OTg)  - G • UN 

where  a small  component  of  gravitation  has  been  dropped  in  the 
expression  for  Lp. 

We  have  to  find  only  one  more  component  of  the  lift  vector, 
namely  Lp,  and  this  is  determined  by  condition  (3) 
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L • V = 0 
If  V is  written 


V = 


Vd 


V u + 
p p 


V u 
N N 


then  LQ  is  seen  to  be 

L V„  + L V 
. _ P P N N 

LD  “ V 

D 

provided  VQ  is  not  zero.  If,  due  to  some  earlier  maneuvering,  VD  is 
zero  or  negative,  the  commanded  lift  vector  is  taken  to  have  the 
direction  Un  and  be  of  some  large,  nominal  value  until  VD  becomes 
positive. 

This  completes  the  guidance  law  derivation.  The  equations  are 
summarized  below  for  convenience.  As  a final  comment,  note  that  a 
multiplicity  of  (false)  target  vectors  and  associated  dive  lines  can 
be  given  sequentially,  if  desired,  for  trajectory  shaping. 


GUIDANCE  EQUATION  SUMMARY 


r *=  Ivj,  - R 
V **  R - a)  x R 


E = r • Up 
H « - V * Up 


F ■ r • 


U 


N 


- V 


r . r 
V • r 
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DISCUSSION  AND  SIMULATION  RESULTS 


In  an  actual  implementation  of  the  guidance  law  derived  in  the 
preceding  section,  a number  of  factors  suddenly  become  of  increased 
importance.  The  most  obvious  one  is  that  when  time  to  go  approaches 
zero,  the  coefficients  A,  B,  C and  D become  singular.  This  is  not 
surprising  since  if  one  has  only  a very  small  correction  to  make  but 
no  time  in  which  to  make  it,  the  required  acceleration  becomes  infi- 
nite. This  difficulty  is  easy  to  overcome  in  the  original  form  of 
the  guidance  law:  simply  stop  recomputing  (freeze)  the  coefficients 
when  time  to  go  becomes  small,  in  the  case  at  hand  somewhere  around 
one  to  one-half  second.  This  is  potentially  important  for  realistic, 
noisy  systems.  It  is  less  clear  how  to  accomplish  this  freezing  in 
the  simplified  form  derived  in  Appendix  B since  indeterminacy  arises 
from  the  cross  product  itself.  In  the  noise-free  cases,  however, 
both  forms  lead  to  impact  errors  which  are  very  close,  differing  at 
most  by  a few  feet. 

The  order  in  which  computations  are  performed  is  important  if  a 
fixed-point  computer  is  used.  Care  can  greatly  reduce  the  dynamic 
range  of  the  variables  involved.  This  must  be  considered  together 
with  the  computation  of  factors  common  to  several  expressions.  If 
airborne  computations  become  a real  problem,  which  seems  rather  un- 
likely with  today's  technology,  it  may  well  be  possible  to  compute 
successive  values  of  the  coefficients  by  perturbation  techniques, 
thereby  reducing  the  computational  load. 

For  the  simulation  results  shown  in  Table  1,  a nominal  reentry 
trajectory  and  three  nonstandard  trajectories  are  considered.  The 
initial  conditions  are  as  shown.  For  the  disturbed  trajectories,  only 
the  condition  which  differs  from  the  nominal  is  given.  Linerized  aero- 
dynamics were  used  for  all  runs.*  The  angular  error  recorded  is  the 
angle  between  the  dive  line  and  the  terminal  velocity  vector.  Differ- 
ences of  only  a few  feet  are  not  considered  significant. 


♦Internal  Memorandum:  Maneuvering  Reentry  Vehicle  Formulation  for 
the  Advanced  System  Simulation,  10  Nov  1977,  Code  CK-21 


INITIAL  CONDITIONS  FOR  SIMULATIONS 


Case  A (Nominal) 

100K  ft  altitude 

200K  ft  up  range  from  the  target 

800  ft  crossrange  to  the  right  of  the  vertical  plane 
containing  the  dive  line 
20K  ft/sec  speed  (in  plane) 

-24°  flight  path  angle 

Lift  vector  in  plane 

Dive  line:  74°  up  from  the  horizontal 

Case  B (Position  Disturbance) 

1800  ft  crossrange 

Case  C (Velocity  Disturbance) 

Velocity  vector  20°  out  of  plane 

Case  D (Acceleration  Disturbance) 

Lift  vector  90°  out  of  plane 
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TABLE  1 


PRELIMINARY  RESULTS 


DOWNRANGE 
MISS  (FEET) 


CROSSRANGE 
MISS  (FEET 


ANGULAR 
ERROR  (DEG.) 


The  miss  distances  shown  in  Table  1 are  trivial  and  demonstrate 
feasibility  of  the  proposed  guidance  equations.  A much  more  detail- 
ed study  considering  noise  effects,  computer  word  length  limitations, 
and  similar  problems  is  required  to  obtain  accuracy  for  a realistic 
system. 


This  appendix  provides  definitions  for  primary  quantities 
involved  in  deriving  the  guidance  equations.  Less  important  terms 
are  defined  as  they  occur  in  the  text. 

COORDINATE  FRAME 


Earth  fixed  dive/ turn: 

Origin  - at  the  target 

X axis  - defined  by  the  unit  vector  UD  parallel  to  the  dive 
line  and  positive  into  the  earth.  The  dive  line  is 
specified  by  the  relative  approach  azimuth  and  flight 
path  angle.  ' 

Z axis  - specified  by  the  unit  vector  in  the  geodetic 
vertical  plane  at  the  target,  positive  into  the 
earth 

Y axis  - specified  by  the  unit  vector  Up  to  complete  right- 
handed  system 

See  Figure  A-l  for  a sketch  of  the  geometry.  The  slight 

distinction  between  geodetic  and  astronomic  vertical  has 

been  ignored  for  simplicity  in  drawing  the  figure. 

Guidance : 

Origin  - center  of  the  earth 

Axes  - inertial,  specified  by  the  guidance  system 


VECTORS 

H - from  the  earth  center  to  the  reentry  vehicle 

Rt  - from  the  earth  center  to  the  earth  fixed  target 

r - line  of  sight  vector  from  reentry  vehicle  to  the  earth 
fixed  target 

V - derivative  of  R with  respect  to  an  earth  fixed  frame 
L - commanded  lift 
G - gravitation  vector 


SCALARS 


E - projection  of  r onto  Up 
H - time  rate  of  change  of  E 
F - projection  of  r onto 
W - time  rate  of  change  of  F 


APPENDIX  B 

CROSS  PRODUCT  FORM  OF  SIMPLIFIED  GUIDANCE  LAW 


Christopher  Gracey 


CROSS  PRODUCT  FORM  OF  SIMPLIFIED  GUIDANCE  LAW 


We  seek  to  show  that  an  equivalent  "cross  product"  form  of  the  guidance 
law  can  be  obtained  when  gravity  is  neglected.  To  that  end  we  start  with  the 
cross  product  form  of  the  guidance  law,  namely 

Ac=  Kajv  x £(F  - K<ji  Tg  V)  x UD]  J B.l 

wi  th 


Ko  = -r-A-— - , K<t  = 2/3. 

ip--u0 


Expanding  the  cross  product,  we  obtain 

Ac  ={(V-UD)(r  - K*  TG  V)  - V*(r  - K TQ  V)UD}  Ka 

Using  the  fact  that 

V = (V-UD)UD  + (V-UN)UN  + (V-Up)Up 

r = (r-U0)UD  + (r-UN)UN  + (r-Up)Up 

We  obtain  for  equation  B.2 


Ac  • KjvuD[(r  - ^ Tg  V)-UDUD  ♦ (r  - K,  Tg  vHyjj, 
Mr  - ^ tg  V)-Vp].  V-(r  - K*  TG  v)u0} 


For  K0  * 2/3,  equation  B.3  becomes 


«C  ■ K» 


V-Wn . V-'Wp 


» V J. 


6/T 


8-1 


B.2 


B.3 


t 


+ V*UD(r  - 2/3  TgV)*UdU0-  V*(r  - 2/3  TgV ) UD 


B.4 


where. 


A,.  = 


-[r-2/3TG  v]  UH 


Ap  * [r  - 2/3  Tg  V 


]'Up. 


This  proves  the  equivalence  of  the  cross  product  form  of  the  guidance  law, 
equation  B.l,  and  the  form  of  the  guidance  law  used  in  the  main  text  when  the 
gains  for  the  cross  product  law  are  properly  selected. 

Rearranging  terms  and  simplifying,  equation  B.4  can  be  written  as 


Af  = K 

L a 


an  v*udun  . Ap  v‘upup  r - - 


6/t; 


6/T 


- [(V-UN)(ruN  - 2/3  T6VUN) 


+ ( V * Up ) ( r-  Up  - 2/3  TGV*Up)]uD 


Ka 


6 fit 


’N'“D“N  ’ "P  ’ UDUP"  [V' 


A„V-U„U„  * An  V-UnU, 


UN  ♦ Ap  VUp]uD 


B.5 


Finally,  equation  B.5  can  be  put  in  the  required  form 

AC  = ^‘UdSANUN  + AP^p  + ADUDj 
' G 


where 


aD  = ' 


apVv  ! Vn’v 

UD*V 


- 

G UD-V 
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